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ABSTRACT 

The damped Lya absorbers (DLAs) in quasar spectra are believed to be the 
progenitors of present-day disk galaxies. We examine the probability for microlensing 
of background quasars by stars in their DLAs. Microlensing by an individual star 
should magnify the continuum but not the broad emission lines of the quasars. 
Consequently, the equivalent width distribution of microlensed quasars would be 
distorted. We model a representative spiral galaxy as a closed system composed of 
a bulge, a disk, and a halo, and evolve the mass fraction of stars in the disk based 
on the observed metallicity of DLAs at high redshifts. The microlensing signatures 
are stronger if the halo of the galaxy is made of Massive Compact Halo Objects 
(MACHOs). In this case, the distortion imprinted by microlensing on the equivalent 
width distribution of quasar emission lines can be detected with high significance in 
a sample of ~ 10 DLAs with HI column densities N ^ 10^^ cm“^ and absorption 
redshifts Zabs ^ 1- About a tenth of all quasars with DLAs {N ^ 10^® cm“^) might 
show excess variability on timescales shorter than hve years. A search for these signals 
would complement microlensing searches in local galaxies and calibrate the MACHO 
mass fraction in galactic halos at high redshifts. 
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1. Introduction 

Damped Lya absorbers (DLAs) are thought to be the progenitors of present-day disk galaxies 
(Wolfe 1995). The abundance of heavy elements at low ionization stages in the absorbers and the 
velocity field traced by them is consistent with values expected for disk galaxies (Turnshek et al. 
1989; Wolfe et al. 1993; Lu et al. 1993; Pettini et al. 1994; Lu & Wolfe 1994). Recent HIRES 
observations on the Keck telescope indicate that the weak, low-ionization, metal absorption lines 
in these systems often show the highest column-density component at the edges of the velocity 
profile, as expected for absorption by a rotating gaseous disk (Wolfe 1995; Prochaska &: Wolfe 
1997). Observations of redshifted 21-cm absorption and emission from DLAs indicate disk-like 
structures of galactic dimensions (Briggs et al. 1989; Wolfe et al. 1992, Briggs et al. 1997), and 
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Faraday-rotation observations are consistent with the existence of micro-Gauss magnetic fields in 
these systems (Wolfe, Lanzetta, &: Oren 1992; Welter, Perry, & Kronberg 1984; Perry, Watson, & 
Kronberg 1993). 

Since the total comoving density of HI in DLAs is comparable to that of stars in the local 
universe, it is only natural to postulate that the cold gas out of which most of the present 
population of stars had formed, was already assembled in galaxies at z ^ 3 (Lanzetta et al. 
1995; Wolfe et al. 1995; Storrie-Lombardi et al. 1996a, 1996b). However, the typical metallicity 
of Z ~ 1O%.^0 (Smith et al. 1996; Lu, Sargent, & Barlow 1996, and references therein) and 
dust-to-gas ratio (Fall & Pei 1993; Pei, Fall & Bechtold 1991) in these systems imply that star 
formation was only at its infancy at these early times. This assertion is consistent with recent 
determinations of the star formation rate at high redshifts (Madau 1996; Lowenthal et al. 1996, 
and references therein). Since star formation requires cold gas and most of the HI detected 
through Lya absorption lines is in damped systems, DLAs are the natural sites for star formation 
at high redshifts (Fall & Pei 1996). In order to unravel the star formation history of the universe, 
it is of fundamental importance to probe the stars and not only the gas in DLAs. 

Direct images of distant DLAs typically reveal a luminous galactic core which is separated 
by 10-20 kpc from the line-of-sight to the quasar (Steidel et al. 1994, 1995, 1996). However, 
any inference about the fraction of gas converted into stars in these systems requires a prior 
knowledge of the Initial Mass Function (IMF) of these stars and their formation history. In 
addition, a major fraction of the baryonic mass might reside in the outer low surface-brightness 
halos of these systems. This possibility is raised by recent microlensing searches (Alcock et al. 
1996), which indicate that a non-negligible fraction of the massive halo of the Milky-Way galaxy 
might exist in the form of Massive Compact Halo Objects (MACHOs). As a supplement to local 
microlensing searches, it would be particularly interesting to examine whether MACHOs populate 
the halos of galaxies at high redshifts. While the microlensing probability is only 10 ® in the 
Milky-Way halo (~ 10 kpc), its value increases up to unity in the cores of halos at cosmological 
distances (~ 5 Gpc). This follows from the linear dependence of the lensing cross-section on the 
observer-lens distance for a source at infinity. Near the center of a DLA, macrolensing by the 
entire galactic potential might take place and yield widely separated quasar images; the likelihood 
of macrolensing and its generic signatures were discussed recently by Bartelmann & Loeb (1996) 
and by Perna, Loeb, & Bartelmann (1997). 

In this paper we quantify the expected microlensing probability in distant spiral galaxies 
which show up as DLAs in quasar spectra. The obvious advantage of these galactic systems is 
that they are selected based on their proximity on the sky to a compact source of light in the form 
of a quasar. Detection of microlensing of the quasar can be used to probe the stellar mass fraction 
of DLAs and to test whether the halos of galaxies at high redshifts are made of MACHOs. The 
characteristic Einstein radius of a solar mass lens at a cosmological distance is ~ 5 x 10^® cm, 
comfortably in between the scales of the continuum-emitting accretion disk (^ 10^^ cm) and the 
broad line region (~ 3 x 10^^ cm) of a bright quasar. This implies that a single microlensing event 
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would magnify the continuum but not the broad lines emitted by the quasar. The lines would 
only be affected by the macrolensing effect of the galaxy as a whole. As a result of this differential 
amplification, the equivalent width distribution of the broad emission lines (Francis 1992) will be 
significantly distorted in a sample of microlensed quasars. Microlensing would also result in excess 
variability of such quasars. In the following, we will quantify the above signatures for a model 
spiral galaxy which acts as a high redshift DLA. 

The existence of a population of normal disk galaxies at high redshifts is supported by recent 
Keck HIRES spectra of the damped Lya absorber towards the quasar Q 2233+1310 (Lu, Sargent 
& Barlow 1997). The metal absorption lines of the absorbing galaxy are shifted relative to its 
Lya emission redshift (Djorgovski et al. 1996), indicating a rotation velocity > 200 km s“^ at a 
separation ~ 20 kpc. Here we limit our attention to galaxies of this population. For concretness, 
we focus our discussion on a single spiral galaxy and extrapolate its present-day properties back 
in time using the simplest closed-box model for its baryonic content. Based on the observed 
metallicity of DBAs, we assume that the stars constitute only 10% of the total disk mass at z ~ 3, 
and that HI is gradually converted into stars until the present epoch. The galaxy is modeled as 
having three components; an old stellar bulge, a disk made of gas and stars, and a halo. We 
compute the significance of microlensing for each of these components. Our simplified model 
ignores the diversity among different disk galaxies or DBAs and focuses on a typical massive spiral 
galaxy. Its only purpose is to examine whether microlensing could be detected in this idealized 
case. Recent imaging of DBA galaxies (Be Brun et al. 1996) indicates that while some systems 
are clearly spirals, others have irregular morphologies. A more diverse treatment of this problem 
would involve many more free parameters, and is beyond the scope of this paper. Our primary 
objective here is to motivate an observational search for these microlensing signatures, rather than 
to explore the entire possible range of parameter space. 

The outline of the paper is as follows. In §2 we describe the distributions of HI and stars in 
our model for a galactic disk. These distributions are then used in §3 to compute the microlensing 
probability as a function of the intercepted HI column density of the disk. In §4 we repeat the 
calculation for a halo made of MACHOs. Sections 5 and 6 relate these results to observable 
signatures. Finally, §7 summarizes our main conclusions. 


2. Distributions of HI and Stars In the Model Disk Galaxy 

The cross-section for gravitational leasing of a high-redshift source (zs > 2) peaks at 
intermediate redshifts, z ~0.2-0.8 (see, e.g. Fig. 1 in Refregier & Boeb 1996). In this redshift 
range, there is compelling observational evidence indicating that disk galaxies have not accreted 
large quantities of matter. The evidence includes the low level of lopsidedness in spiral galaxies 
(Zaritzky & Rix 1996), the thinness and coldness of the stellar distribution in galactic disks (Toth 
& Ostriker 1992), the small scatter in the Tully-Fisher relation (Eisenstein & Boeb 1996), the 
similarity between the kinematics at a given luminosity of galaxies out to z ~ 1 and local galaxies 
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(Vogt et al. 1996), and the statistics of multiply imaged quasars (Mao &: Kochanek 1994; Rix et 
al. 1994). We therefore anchor our modeling to the current properties of a typical galactic disk, 
and evolve the gaseous mass fraction of its disk based on a simple closed box model. 

Observations of local spiral galaxies show that their HI gas is generically distributed in an 
exponential disk with a radial face-on column-density profile (Broeils & Von Waerden 1994), 

^face-on(^^ now) = Nq exp(-R/Rg), (1) 

where Nq is the central column density and Rg is the scale radius of the HI disk. For our model 
galaxy we adopt Nq ~ 10^^cm“^, which is the typical value observed in nearby galaxies. 

The surface mass density of stars in a nearby spiral galaxy can be modeled as the sum of an 
exponential disk and a bulge-described by a de Vaucouleur law, 

5 .face-on(^^ now) = Sq exp(-R/Rs) + Sg exp{-7.67[(R/Re)3 - 1]}, (2) 

where R is the radial coordinate on the face of the disk. The central mass density of stars in the 
disk, Eq, can be calibrated from the characteristic surface brightness observed in galactic disks 
/q ~ 21.65^s, where hb is B mag arcsec“^ (Freeman 1970). We convert the surface brightness 
to mass density assuming a mass to light ratio of 5 in solar units (Broeils &: Courteau 1996). 
The bulge surface density at 72 = i?e, Eg, can be related to Eq and the total disk and bulge 
luminosities, Lj and Lb, through the relation Ld/Lb = 0.28(72s/72g)^(Eo/Ee). This relation results 
from an integral of the corresponding surface brightnesses over the area of the disk. For typical 
values of Ld/Lb ~ 1.3 and Re/Rs ~ 0.4 (Freeman 1970; Burstein 1979), we get Eg/Eg ~ 1.5. The 
scale radius of the stars Rs is strongly correlated with the scale radius of the gas Rg. Figure 8b 
in Broeils & van Woerden (1994) illustrates that the HI surface densities of different galaxies fall 
off with very similar profiles when the radius is expressed in units of the optical radius 7225 of 
the 25 mag arcsec“^ isophote in the B band. With 7o ~ 21.65/iB, we derive R 25 ~ 3.172s and 
from the slopes of the profiles in Broeils & van Woerden (1994) we estimate 72g « 0.87225 so that 
72s ~ 0.472g. For our model galaxy we adopt the typical value of 72g = lOkpc. Our results for 
microlensing by stars in the disk are not very sensitive to the particular choice of 72g and mostly 
depend on the adopted ratio between 72g and 72s, for which the observed scatter is small. 

We evolve the model galaxy back in time as a closed system. We assume that the model disk 
was already assembled at z ~ 3 and that its HI content evolved only as a result of star formation; 
this assumption is consistent with the fact that the comoving HI density of DBAs at z ~ 3 is 
comparable to the mass density of stars in the local universe. As nearby galaxies exhibit a much 
older stellar population in their bulges than in their disks, we keep the bulge properties fixed in 
time and evolve only the stellar mass fraction of the disk with redshift for z 3. Based on the 
observed metallicity of DBAs (Pettini et al. 1994, Pettini et al. 1996) we assume that at a redshift 
2 ; ~ 3 the mass fraction of stars in the disk was about a tenth of its present value. Under these 
assumptions, equation (|I|) is generalized to 

+ ANo{z)e-^/^, 


( 3 ) 
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where in the closed box model the HI and stellar surface densities obey 

^.face-on {R,z)= COnst (4) 

with Shi = and therefore 

AShi(-R, z) = S*^disk(-R, 0) - S*^disk(-R, z). (5) 

Here mn is the mass of a hydrogen atom, and the subscript {*, disk} refers to the disk stars only. 
In the context of the closed box model, the gas out of which the disk stars formed must have 
had the same radial distribution as those stars have today, and so equation (^) can be regarded 
as a condition on the corresponding central densities. We constrain the redshift evolution of Shi 
by the condition S* disk(-2 = 3) ss 10% x S* disk(2 = 0), and by the observed redshift evolution 
of the HI content of DLAs (Wolfe et al. 1995; Storrie-Lombardi et al. 1996a). Both of these 
constraints are reasonably satisfied in the redshift interval 0 ^ z ^ 3, if one uses a power-law 
redshift evolution of the form, AShi(I?, z) = m}iNz^ ex.p{—R/R^), with 5 ~ 6 and N = lO^^crn”^. 
Based on equation (^), we then derive 

^*,disk{R,z) = Po - mHAAo(z)] exp(-ii/i?s) = S}(z) exp(-i?/i2s), (6) 

and thus the stellar surface density profile is given by 

S^e-°n(i?,^) = sUz)exp(-i?/i?s) + Seexp(-7.67[(i?/i?e)3 - 1]). (7) 


Equations (^) and (^) describe the radial profiles for a face-on disk. If the normal to the 
disk is inclined by an angle 9 with respect to the line of sight, then the observed distributions are 
changed to 


N{R,e,z) 


[Noexp{-R/Rg) + AjVo(z)exp(-i?/i;s)] 
cos 6 


( 8 ) 


and 


U^ (pj Z) — ^★,disk ' '^★,bulge — ^ 




(9) 


where b = Ryl — sin^ 0cos^ cj) is the impact parameter and (j) is the azimuthal coordinate on the 
sky of the line-of-sight relative to the center of the galaxy. In this derivation we have assumed, for 
simplicity, a razor-thin disk and a spherical bulge. 


3. Microlensing by Stars in the Disk and the Bnlge 

For a single galaxy, the optical depth to microlensing is given by the dimensionless surface 
mass density of its stars 


-'crit 


K* = 


( 10 ) 
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as normalized by the critical surface density, 


^crit — 




dvrGDd-Dds 


( 11 ) 


Here G is the gravitational constant, c is the speed of light, and Dg and Hds are the angular 
diameter distances between the observer and the lens, the observer and the source and lens and 
the source, respectively (Schneider, Ehlers, & Falco 1992). Throughout the paper we assume 
Hq = 50 km s“^ Mpc“^, Hq = 1; and a quasar redshift of Zg = 4. 


We would like to find the probability distribution for the microlensing optical depth, P(fi;*), 
given a particular observed value for the HI column density in our model galaxy. We assume 
random inclinations of the galactic disk relative to the line-of-sight. Given the HI column density 
in equation (|^, the inclination angle is restricted by the condition 9 > 9q, where 


cos 6*0 = 7 o = min 



No + ANo{z) \ 
N ) • 


( 12 ) 


The uniform probability distribution for the cosine of the inclination angle is then 

P(cos 9)d cos 9 = — d cos 9 for 0 < cos 0 < 7 o- 
70 


(13) 


Given an inclination angle, 9, the microlensing optical depth k* can be obtained from 
equations dH) and (jl^). The value of k* depends on the observed column density N through the 
radius R intercepted by the line-of-sight in the plane of the disk; this radius is in turn obtained by 
inverting equation (^). Given an observed HI column density iV, the probability distribution for 
the optical depth k* is 


P{^^AN) = 


1 


r2-n 


27r7o Jo 
1 


dcos 9 


r2TT 


with 


27 r 7 o Jo 

dR 


dK-t 


-imbK^h-^Rl 


dh 


dcos 9 
N 


Kd 

cos'^9 


1 + 


cos 9 dR 
Rg d cos 9 


-1 


(14) 


(15) 


dcos 9 Rg ^No ex.p{—R/Rg) -|- Rg ^ANo{z) ex.p{—R/Rs) ’ 

where ko = Sbuige/^crit and Kd = Sdisk/Scrit- Note that P(k*|A*') depends on the absorber redshift 
through equations (0) -(|I|) and (H). 


If in total K ^ 1, then the galaxy as a whole produces multiple images. Under these 
circumstances there could still be microlensing, as observed in the lens Q2237-I-0305 (Irwin et 
al. 1989; Wambsganss, Paczyiiski, &: Schneider 1990; Rauch &: Blandford 1991; Racine 1991). 
However, in this case the observed spectrum will not show a single column density since different 
quasar images intercept the disk at different impact parameters and therefore probe different HI 
column densities. The Lya absorption feature will then be composed of multiple troughs, each 
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having a depth proportional to the magnification of its corresponding image and a width scaling 
as the square root of the HI column density probed by this image (Bartelmann & Loeb 1996; 
Loeb 1997). Since the interpretation of N is not unique in this case and since most DLAs with 
N ^ 10^^ cm“^ are not macrolensed (Bartelmann & Loeb 1996), we restrict our attention to 
values of K* ^ 1. Figure shows the probability at a given value of N for obtaining a value of k* 
smaller than unity but higher than some threshold Kmim < k* < 1 |A^) = . dK*P(K*|A'). 

The different panels correspond to different redshifts of our model galaxy; panel (a) refers 
to Zabs = 0.5 and panel (b) corresponds to Zabs = 2.5. For a low absorber redshift, the 
probability of having a microlensing optical depth k* > Kmin = 0.05 is negligible at low column 
densities {N 3 x 10^^cm“^), but becomes significant above 10^^cm“^. For a high absorber 
redshift, the probability of having k* > 0.05 becomes significant at yet higher column densities, 
N ;^4-9xl0^^cm“^. This is due to the increase in Sent as the lens gets closer to the source. In 
order to reach the same values of k* at high redshifts, it is necessary that the line-of-sight to the 
quasar will pass closer to the center of the absorber, where the HI density is higher. The drop 
in the probabilities at yet higher HI densities is due to the increase in the probability of finding 
systems with k> 1 (which we exclude from our analysis). 


4. Microlensing by MACHOs 

We now consider the possibility that in addition to the stars in the disk and the bulge, there 
is a population of compact objects which dominate the mass of the halo in our model galaxy. 

We model the spherically-averaged mass distribution of the entire galaxy as a singular 
isothermal sphere (SIS) with a one-dimensional velocity dispersion cr. 


, . 2cr^ 
Msis(r-) = —r, 


(16) 


where r is the radial spherical coordinate and a the velocity dispersion for which we adopt a 
typical value of 170 km/sec. The SIS model yields a flat rotation curve, as observed in nearby 
spirals. We assume that the total mass of the galaxy is composed of the disk of gas and stars [cf. 
Eqs. and d^] and a spherical component (denoted by SP) which includes both the bulge and 
the halo, 

Msis{r) = Msp(r) + MHi,disk(?’) + M*,disk(?’)- (17) 

The HI and stellar masses out to a radius R = r on the face of the disk are found by integrating 
their respective surface densities, 


A7HI,disk(?’) = 


2tt / 


HI 


ex.p{—R/Rg)RdR 




1 - exp(-r/i?g) 1 + 


R^ 


(18) 




with Syj = TTiH-ZVo and 




27rSo-Rs^ 


1 - exjp{-r/Rs) 



(19) 


The mass density of the spherical component (bulge plus halo) can be derived from equation (0) 
through the relation, 


/ 0 sp(?’) = 


1 d-Msp 
dvrr^ dr 


a 


2TTGr‘^ 


1 

Jr L 


SHiexp(-r/it’g) + Soexp(-r/ij’s) 


( 20 ) 


The surface density as a function of the impact parameter on the sky b is then found by integrating 
this mass density along the line of sight S( 6 ) = d(p{b, () with r = ( 6 ^ + C^) 2 . This yields 


Ssp( 6 ) 


2Gb 


r(i/2) 


Ko 



+ So Ko 



( 21 ) 


where Kq is the modified Bessel function of the zeroth order. The optical depth to microlensing is 
then given by 


k{R, 9, (j)) 


KSP + «^*,disk 
1 


-'crit 


Ssp(fe) + 


Soexp(-i?/i?s) 

cos 9 


( 22 ) 


where, as before, 9 is the inclination of the normal to the disk relative to the line-of-sight, 4> is 
the azimuthal coordinate on the sky, R is the radial coordinate in the plane of the disk, and 
b = R^J 1 — sin^ 9 cos^ </>. The probability distribution of k* for an observed HI column density N, 
can be calculated again as in equation (jlj), 


P(k*|A^) = 


db 


1 


r2Tr 


27r7o Jo 


dcos 9 


SIS 


62 


R. 




R(r 


^*,disk 

COs 20 


1 + 


COS 9 dR 
d cos 9 


-1 


(23) 


where Ki is the modified Bessel function of the first order, and dR/dcos 9 is given in equation (IT). 


Figure ^ shows the probability for obtaining a value of k* smaller than unity but higher 
than some threshold Kmin at either a low absorber redshift [z^hs = 0.5 in panel (a)] or a high 
absorber redshift [z^bs = 2.5 in panel (b)]. As evident from comparing Figures 2 and 1, the 
existence of MACHOs enhances considerably the microlensing probability in DBAs which are not 
macro-lensed. This is especially true for absorbers with lower column densities. 


5. Observational Signatures of Microlensing 

The microlensing probability derived in the previous sections can be related to two observable 
signatures: (i) distortion in the equivalent-width distribution of the broad emission lines of 
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microlensed quasars (Canizares 1982), and (ii) excess variability of microlensed quasars (see, e.g. 
Wambsganss &: Kundic 1995). In the following we quantify both of these effects as a function of 
the observed HI column-density for our model galaxy. 


5.1. Equivalent—Width Distribution of Microlensed Quasars 

Quasars can be significantly affected by microlensing only if the size of their emission region 
is smaller than the projected size te of the lensing zone, i.e. the “Einstein radius” of the lensing 
star. The maximum magnification of a finite circular source of radius and uniform surface 
brightness is given by [Schneider et al. (1992)] 

Itmax = Y^l + 4(rE/rs)2 . (24) 

At cosmological distances, the lensing zone of a star of mass bas a characteristic scale 

of ~ 5 X 10^® (Mstar/.^©)^^^ cm. In comparison, the optical continuum emission in quasars is 
believed to originate from a compact accretion disk. The UV bump observed in quasar spectra 
is often associated with thermal emission from an accretion disk with a surface temperature 
Tdisk = 10® Ts K, where Ts 1 (e.g. Laor 1990), and so the scale of the disk emission region 
must be ~ 10 ^® cm , where L 46 is the corresponding luminosity of the quasar in units 

of 10^® erg s“^. Thus, for lens masses Mg^ar S> 1 O“®M 0 , the continuum source is much smaller 
than the lensing zone and could therefore be magnified considerably. This expectation is indeed 
confirmed in the nearby lens of Q2237+0305, where variability due to microlensing has been 
observed (Wambsganss et al. 1990; Rauch & Blandford 1991; Racine 1991; see also Gould &: 
Miralda-Escude 1996). However, reverberation studies of the time lag between variations in the 
continuum and the line emission in active galactic nuclei indicate that the broad emission lines 
of quasars originate at a distance of ~ 3 X 10^^ cm (Peterson 1993; Maoz 1996). For a 
solar mass lens with luminosity 10^® erg cm“^, use of Equation ( p4|) shows that the maximum 
magnification differs from unity by only a few percent. This implies that microlensing of the 
broad line region by a single star can be safely ignored. The broad line region would only be 
macrolensed by a factor (/r) due to the average effect of all the stars. As a result of the difference 
in amplification between the lines and the continuum, the equivalent width of the lines will change 
during a microlensing event. For compact objects which are uniformly distributed throughout the 
universe, the microlensing optical depth is a function of redshift and should therefore result in an 
apparent evolution of the equivalent-width distribution of quasars. The lack of such evolution in 
existing quasar samples was therefore used to set an upper limit on the mean cosmological density 
of stars and MACHOs in the universe (Dalcanton et al. 1994). Nemiroff (1988) examined the 
effect of microlensing on the shape of the broad emission lines, and Hawkins (1996) addressed the 
potential to detect cosmologically distributed microlenses based on quasar variability data. Here 
we propose to look for microlensing signatures through a comparison between the equivalent width 
distribution of quasars with DLAs and those without. Since quasars with DLAs are selected based 
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on their small projected separation from foreground galaxies, they are more likely to be lensed by 
stars or MACHOs than the rest of the quasar population. The latter approach is likely to have a 
larger signal-to-noise ratio if most of the stellar objects in the universe are grouped into galaxies. 
Since DLAs appear only for a small fraction of all quasars, ~ 10%(A^/10^‘^ cm it should be 

straightforward to obtain an accurate equivalent-width (EW) distribution for a control sample of 
quasars which are not microlensed. 


In order to find the effect of microlensing on the EW distribution of quasar emission lines, 
it is necessary to calculate the magnification probability as a function of k*. This probability 
function was derived in the literature through extensive numerical simulations (Witt 1993; 
Wambsganss 1990; Lewis et al. 1993; Kundic & Wambsganss 1993; Wambsganss &: Kundic 
1995). In particular, the simulations of Rauch et al. (1992) demonstrated that the probability 
distribution of magnifications for a point source seen through a field of stars can be reasonably 
well approximated by the analytic expression proposed by Peacock (1986) 

= (25) 

a/r 

where 

= 2x[l - 1^(1^^ - 1)"^'^^] and df/d^ = 2x(/i2 - l)-^/^ 
and the parameter x is related to the dimensionless surface mass density k* through the equation 

roo I 

{n)= /rP(Ai,x)dM = 7—-(26) 


The probability distribution of magnifications for a given value of the observed HI column density 


N can be obtained by combining the probability distribution for k* at a given N [cf. Eqs. (IT) or 
with the probability distribution of magnifications at a given k* [cf. Eq. 


Pifi\N)= [' dK^P{^\K,)P{K,\N). 

JO 


(27) 


We define /j/(Ao) to be the intensity of the quasar continuum in the neighborhood of the 
wavelength Aq of a particular emission line and AI^(X) to be the difference between the total 
measured intensity and the continuum. If the continuum is magnified by a factor //, then its 
observed intensity changes to ^Iy{Xo). The extended broad line region is affected by the combined 
effect of many stars, and so the intensity of the lines is enhanced by the average magnification 
factor (^) in equation (^). Consequently, the equivalent width (EW) of the emission line, defined 
as 




Ah 

Iu{^o) 


dX 


is changed by a factor (/r)/^, namely 


fPA = VEo— 


(28) 


(29) 
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where Wq is the intrinsic EW of the unlensed quasar. 

Even in the absence of lensing, quasars do not possess a unique EW value in their emission 
lines but rather show a wide probability distribution of EW values (Francis 1992), which we define 
as P{Wo). The magnification due to lensing distorts this probability distribution. The modified 
distribution of EW values for quasars which show DLA absorption with a column density N would 
be, 

PiWx\N) = d^xP{^^\N)P (30) 

The observed probability distribution of EW can be well approximated as a log-Gaussian function 
(see also Dalcanton et al. 1994), 

P(lTo) = exp{-[ln(lTo) - uj]V2a^w}, (31) 

^j2na‘^Wo 

where the parameters uj and aw obtain different values for different emission lines. We optimize 
these values so as to fit the data from the Large Bright Quasar Survey (LBQS) reported by Francis 
(1992). 

Figures 3 and 4 illustrate the unlensed distributions (solid lines) chosen to best represent 
the observational data from LBQS, part of which is shown as squares. The figures also show 
the distorted EW distribution for quasars which are lensed by stars in the disk (dotted-dashed 
line), disk+bulge (dashed line) and disk+bulge+halo (dotted line). Figure 3 corresponds to 
our model galaxy being at a low redshift, Zabs = 0-5 and with an observed column density 
N = 5 X lO^'^cm"^ (Fig. 3a) or = lO^^cm”^ (Fig. 3b). Figure 4 shows the results for z^hs = 2.5 
with N = 5 X lO^^crn”^ (Fig. 4a) or A" = 9 x lO^^cm”^ (Fig. 4b). Note that the area below 
the curves is in some cases smaller than unity, because we restrict our attention to systems with 
/t* < 1. In all cases, the distortion of the EW distribution is more pronounced for galaxies with 
MACHOs. The level of this distortion can therefore be used as a sensitive probe of the composition 
of galactic halos at high redshifts, and by that complement searches for MACHOs in the local 
universe (Alcock et al. 1996, and references therein). 

In order to estimate the minimum size of the quasar sample that allows detection of the 
distortion to the EW distribution, we performed a numerical simulation of random deviates drawn 
from the lensed (i.e. distorted) distribution and computed the statistic with respect to the 
observed data for the unlensed distribution (Francis 1992). We assumed that the number of 
quasars with DLAs is much smaller than the number of quasars without DLAs, which were used 
to define the unlensed distribution. In this case, ~ ni)‘^/ni], where hi is the number 

of simulated events from the lensed distribution in the i-th bin and rii is the number of expected 
events according to the unlensed distribution. For low-redshift (^^abs ~ 0.5) absorbers with halos 
made of MACHOs and observed column densities N ~ lO^^crn”^, we find that a relatively small 
sample of ~ 10 systems drawn from the lensed distribution can be rejected as being drawn from 
the unlensed distribution with a probability of 5 x 10“^. 
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If MACHOs are absent, the number of systems required to get the same confidence level is 
increased to ~ 20. For high redshift (zabs ~ 2.5) absorbers with MACHOs and N ~ 10^^cm“^, 
the required number of systems, ~ 100, exceeds the number of all known DLAs. However, for 
high-redshift absorbers with column densities on the order of ~ 8 — 10 x 10 ^^ cm“^, even a small 
sample of ~ 10 systems could be sufficient to detect the signal at the same confidence level. 

We note that our calculations were formulated in terms of the stellar surface density k*, 
and therefore ignored contributions to the magnification distribution from shear 7 , or additional 
sources of k (e.g. due to a smoothly distributed dark matter component). The associated increase 
in the average magnification (^) would affect both the lines and the continuum. In addition, the 
finite source size would tend to reduce the expected signal relative to our point source calculation. 
In reality, the magnification is limited to a value fimax ~ for <C (see Equation (^). For 

a stellar mass of ~ O.SMq and a continuum emitting region of a quasar of size ~ 10^^ cm, we find 
IJ'max ~ 10^. The contribution to F’(//) from higher values of // is negligible; P{^ ^ 100) ^ 10“^ 
for K* = 0.1 and P{^ ^ 100) ^ 10“^ for k* = 0.8. 

Our calculation of the EW distribution ignores any magnification bias due to lensing, by 
which faint QSOs are raised above the detection limit. This bias could affect our analysis only if 
the EW were to depend on the luminosity of the source. Such a dependence, otherwise known as 
the “Baldwin Effect” (Baldwin 1977; Baldwin et al. 1988), is generally very weak and was only 
marginally detected for the CIV line (Francis 1992). By modeling EWoc L~^'^ , this weak scaling 
results in a ~ 1% change of the EW values (due to magnification bias) for k* = 0.1 and in a ~ 5% 
change for k* = 0.8. Thus, the magnification bias can be safely ignored in our calculation. 


5.2. Time variability 

The microlensing signature discussed in the previous subsection can be detected at one point 
in time, as it only requires a statistical comparison between the EW distribution of quasars which 
have DLAs in their spectra and that of quasars which do not. By monitoring a quasar over a 
sufficiently long period of time, it is possible to search for yet another signature of microlensing, 
namely an excess temporal variability of the quasar flux due to the motion of the microlenses 
relative to the quasar. Below we quantify the probability for detecting this excess variability as a 
function of the observed HI column density in our model galaxy. 

The net motion of the stars relative to the quasar is the sum of their motion within the 
galaxy and the bulk velocity of the galaxy relative to the line-of-sight to the quasar. Wambsganss 
& Kundic (1995) have calculated the probability distribution for observing an event with duration 
smaller than the crossing time of the stellar Einstein-radius to as a function of k* for a population 
of single-mass stars. Given the source and lens redshifts, the value of to depends on the mass 
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Mstar of the lensing star and on its transverse velocity v relative to the line-of-sight, 


to — 


^{4GM,t^r/c^){DdDds/D, 


(32) 


while an “event” was defined as an increase in the magnification by a factor Am ^ 0.2 mag. 
Wambsganss &: Kundic (1995) presented the probability distribution of to for simulations with 
K* = 0.2 in their paper, and have kindly provided us with their results also for k* = 0.5 and 
K* = 0.8 (Wambsganss &: Kundic 1996, private communication). In order to obtain a rough 
estimate for the probability distribution P{t < tmaxl^*); we have interpolated these numerical 
results over the entire range of 0 < k* < 1. Then, by combining this distribution with equation (14) 


or equation (23), we have obtained the probability for observing a microlensing event with duration 
less than tmax given that the observed HI column density is N, 


P{t < tn 


C \N) = ['dK^P{t< 

^max \k,)P{k,\N), 

Jo 


(33) 


The transverse velocities of stars in galaxies are typically of order 150-200 km s“^. The 
characteristic transverse velocity of galaxies due to large-scale structure can be as large as 
500 km/sec (Strauss &: Willick 1995, and references therein). We therefore parametrize, 

V = 300 X usoo km/sec. The simulation of Wambsganss & Kundic (1995) assumes a single mass for 
all stars. We estimate the numerical value of to by averaging over a Scalo (1986) mass function for 
the stars, with (M^/^) « O.63M0^/^. Using a source redshift Zg = 4, we then get to ~ 28 u3qq yr for 
a lens at Zabs ~ 0.5 and to ~ yr at Zabs ~ 2.5. Note that the above time scales are averaged 

over the stellar mass function. In reality, the low mass tail of the stellar distribution could give 
rise to events of much shorter duration. 

Our results are shown in Figures ^ (for disk+bulge) and ^ (disk+bulge+halo) for a low 
redshift absorber [panel (a)] or a high redshift absorber [panel(b)]. The figures show the probability 
of observing a microlensing event with duration less than t ma x = 0.08to (solid line), 0.12to (dotted 
line), 0.16to (dashed line), and 0.2fo (dotted-dashed line). 

The variable magnification during a microlensing event could be detected through two effects: 


1. Changes in the EW of the broad emission lines with time. A magnification factor of Am 0.2 
mag will reduce the EW by a factor 0.2 over time. This change could be detected through 
a monitoring program of quasar emission lines, of the type realized in reverberation studies 
of quasar variability (Kaspi et ah 1996). It should be generally distinguishable from changes 
due to intrinsic variability, where an increase in the continuum intensity is followed, after a 
certain time, by a corresponding increase in the intensity of the emission lines (Maoz 1996). 

2. Variability in the QSO light curve. Excess variability in the quasar luminosity due to 
microlensing could in principle be distinguished from intrinsic quasar variability. In the low 
optical depth regime (k <C 1 ), the microlenses are isolated and the light curves they produce 
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have generic shapes which are symmetric in time (Paczyhski 1986). Even at higher optical 
depths, there should still be a statistical symmetry between the rising and the falling parts 
of the light curves. Finally, as long as the continuum emission region is much smaller than 
the size of the stellar leasing zone, the excess temporal variation of the quasar should be 
achromatic. 


6. Conclusions 

We have shown that the subset of all quasars which have damped Lya absorption in their 
spectra due to an intervening spiral galaxy should possess an enhanched tail in the equivalent 
width distribution of their broad emission lines, and exhibit excess variability, relative to the rest 
of the quasar population. Both of these microlensing signatures are more pronounced if the halos 
of DBAs are composed of compact objects. 

Previous attempts to find a redshift evolution in the equivalent-width distribution of a large 
sample of quasars due to microlensing by a population of intergalactic stars have yielded a negative 
result (Dalcanton et ah 1994). In this paper we have shown that the significance of microlensing 
relative to the statistical noise should be much more pronounced in a subset of all quasars which 
are located behind galactic HI disks. Using a simple closed-box model for a spiral galaxy we have 
found that disk+bulge microlensing could be detected through its imprint on the equivalent-width 
distortion with a signal-to-noise ratio ^ 2 in a sample of ~ 20 damped Lya absorbers (DLAs) at 
^abs ^ 1 with N 10^^ cm The necessary sample size is reduced to ~ 10 if the galaxy halo 
is composed of MACHOs with a velocity dispersion of a ~ 170 km/sec. The necessary sample 
size should scale as oc for other galactic systems and therefore massive galaxies are likely to 
dominate the statistics. In addition, we find that about a tenth of all quasars with DLAs are likely 
to show excess variability on timescales shorter than hve years (cf. Figures 5 and 6). 

Unfortunately, the current sample of ~ 80 quasars with DLAs (Wolfe et al. 1995) includes 
only several absorbers with z ^ 1 and might not be sufficiently large to demonstrate the existence 
of microlensing. The number of known DLAs could increase by an order of magnitude as a result 
of spectroscopic follow-ups on the catalog of ~ 10^ quasars which is expected to be compiled by 
the forthcoming Sloan Digital Sky Survey (Gunn &: Knapp 1993). Even before any absorption 
data is reduced, it would be interesting to select background quasars which are projected close to 
foreground galaxies and therefore are likely to be microlensed. Related studies by Webster et al. 
(1988) showed an increased number of quasar-galaxy pairs than expected from random alignments. 
The enhancement in the quasar surface density near galaxies was interpreted as macro-lensing and 
used to draw conclusions about the distribution of matter around galaxies. 

Detection of the level of microlensing for quasars with damped Lya absorption can be used 
to calibrate the mass fraction in the form of massive compact objects in galactic halos at high 
redshifts. Figures 3-6 show that the microlensing signal is enhanced when an isothermal halo made 
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of MACHOs is added to the galactic disk and bulge. Direct imaging of the DLAs can be used to 
infer the projected separation between the luminous center of the absorber and the quasar. In 
cases where the DLA redshift is known (e.g. through its Lya emission; see Djorgovski et al. 1996, 
and Lu et al. 1997), it might also be possible to infer spectroscopically the velocity dispersion of 
the intervening galaxy. When combined with the information gathered by microlensing searches 
in the local universe, such studies could extend our knowledge of the composition of galactic halos 
out to redshifts as high as z ~ 5. 

We thank an anonymous referee for insightful comments that improved the presentation. This 
work was supported in part by the NASA ATP grant NAG5-3085 and the Harvard Milton fund 
(for AL) and by a fellowship from the university of Salerno, Italy (for RP). 
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Fig. 1.— Probability of finding the microlensing optical depth k* in the range between ^min and 
1 as a function of the observed HI column density in our model galaxy. Results are shown for an 
absorber (lens) redshift Zabs = 0.5 in panel (a) and Zabs = 2.5 in panel (b). In both panels we show 
three curves for Kmin = 0.05 (solid line), Kmin = 0.1 (dotted line), and ^min = 0.5 (dashed line). 
The source redshift is fixed at Zg = 4. 
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Fig. 2.— Same as Figure 1, but including a halo of MACHOs in addition to the bulge and the disk 
stars. The total microlensing optical depth is k* = /i*,disk + [cf- Eqs. (22) and (23)]. 
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Fig. 3.a.— Effect of microlensing on the equivalent-width distribution of quasars with DLAs. The 
solid line shows the analytical log-Gaussian fit to the observational data (squares) based on the 
Large Bright QSO Survey (Francis 1992). The other lines show the distortion of this distribution 
due to microlensing by stars in the disk (dashed-dotted line), the disk-fbulge (dashed line), and the 
disk+bulge-fhalo (dotted line) of our model galaxy. The galaxy is assumed to be observed with an 
HI column density N = 5 x 10^®cm“^ and a redshift Zabs = 0.5. The source redshift is Zg = 4. 
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Fig. 3.b.— Same as Figure 3a, but with N = lO^^cm 



- 22 - 


Fig. 4.a.— Same as Figure 3a, but with Zabs = 2.5 and N = 


X lO^^cm 
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Fig. 4.b.— Same as Figure 4a, but with N = 9 x lO^^cm 
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Fig. 5.— Probability for observing a magnification of 0.2 mag with duration shorter than tmax as a 
function of the observed HI column density. Our model galaxy is located at a redshift Zabs = 0-5 in 
panel (a) and Zabs = 2.5 in panel (b). In both panels, tmax = O.OSto (solid line), tmax = 0.12to (dotted 
line), tmax = O.lOto (short dashed line), tmax = 0.2to (long dashed line). Here, to ^ 28u^qq yr at 
Zahs = 0.5, and to ^ 17i'^oo ^abs = 2.5, where U 300 is the transverse velocity of the microlenses 
relative to the source in units of 300 km s“^. These numbers were obtained for the average stellar 
mass in a Scalo (1986) mass function. The source redshift is fixed at Zg = 4. 
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Fig. 6.— Same as Figure 5, but including a halo made of MACHOs. 



